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Bulk-boundary correspondence in (3-|-l)-dimensional topological phases 
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We discuss (2+l)-dimensional gapless surface theories of bulk (3+l)-dimensional topological 
phases, such as the BF theory at level K, and its generalization. In particular, we put these theories 
on a flat (2+1) dimensional torus parameterized by its modular parameters, and compute the 
partition functions obeying various twisted boundary conditions. We show the partition functions 
are transformed into each other under SL{3, Z) modular transformations, and furthermore establish 
the bulk-boundary correspondence in (3+1) dimensions by matching the modular S and T matrices 
computed from the boundary field theories with those computed in the bulk. We also propose the 
three-loop braiding statistics can be studied by constructing the modular S and T matrices from 
an appropriate boundary field theory. 
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The bulk-boundary correspondence is one of the most 
salient features of topologically ordered phases of mat¬ 
ter. In topologically ordered states in (2+1) dimensions 
[(2+l)d], all essential topological properties in their bulk 
can be derived and understood from their edge theo¬ 
ries, such as quantized transport properties, properties of 
bulk quasiparticles (fractional charge and braiding statis¬ 
tics thereof), and the topological entanglement entropy, 
etc.^^® Edge or surface theories also play an important 
role in symmetry-protected and symmetry-enriched topo¬ 
logical phases. 

The purpose of this paper is to study the bulk¬ 
boundary correspondence in the simplest (3+l)d topo¬ 
logical field theory, the BF topological field theory^®^^^, 
and its generalizations. The BF theory describes, among 
others, the long wave-length limit of BCS superconduc¬ 
tors, and the deconfined phase of the Zk gauge theory. 
It is also relevant to the hydrodynamic description of 
(3+l)d symmetry-protected topological (SPT) phases in¬ 
cluding topological insulators and related systems. 

To put our purpose in the proper context, let us give 
a brief overview of the bulk-boundary correspondence in 
(2+l)d topologically ordered phases. For (2+l)d topo¬ 
logical phases, bulk topological phases can be character¬ 
ized by the modular S and T matrices. The S and T 
transformations generate the basis transformation in the 
space of degenerate ground states, which appear when 
the system is put on a spatial two-dimensional torus. 
Combined together, the S and T transformations form 
the group S'L(2,Z), the mapping class group of the two- 
dimensional torus T^. Their geometric meanings are the 
90° rotation and Dehn twist defined on the torus, respec¬ 
tively. In the basis in which the T matrix is diagonal (the 
so-called quasi-particle basis), the diagonal entries of the 

matrix encode the information on the topological spin 
of quasi-particles. On the other hand, the S matrix con¬ 
tains the information of the braiding and fusion. For an 
Abelian topological phase, the elements of the S matrix 
are given by braiding phases between quasiparticles, up 
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to an over all normalization factor l/I?, where T) is the 
total quantum dimension. 

On the other hand, at their boundary (edge), gapless 
boundary excitations supported by a (2+l)d topologi¬ 
cal phase can be described by a (l-l-l)d conformal field 
theory (CFT).^^ There is one-to-one correspondence be¬ 
tween quasi-particle excitations in the bulk and primary 
fields living on the edge. On the (l-l-l)d spacetime torus, 
one can form the character Xji'^) from the tower of states 
built upon a primary field Oji 

X,{T)=Trn, (i) 

where Hq and Pq are the Hamiltonian and the mo¬ 
mentum operators, respectively, the complex parameter 
T = Ti iT 2 is the modular parameter parameterizing 
the spacetime torus, and the trace is taken over all states 
in the Hilbert space Hj that is built upon the highest 
weight state associated with the primary field Oj. The 
characters Xj transform into each other under the mod¬ 
ular transformations of the spacetime torus. Under the 
modular T and S transformations, the characters Xjij) 
transform as 

X,(r + l) = e^"('‘^-«)x,(r) = 7-,x,(r), 

= '^SjpXi'iT), ( 2 ) 
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where the matrices T and S represent the action of the 
T and S modular transformations on the characters, re¬ 
spectively. The matrix T is a diagonal matrix and in¬ 
cludes the conformal dimension h for each character and 
the central charge c for the CFT. The T and S matri¬ 
ces for the characters in the edge theory have the direct 
correspondence (and are essentially identical) to the the 
T and S matrices defined for the corresponding (2 -|- l)d 
bulk topological theory. 

Coming back to our main focus, i.e., (3 -I- l)d topolog¬ 
ically ordered phases, the bulk topological system can 
be defined on a spatial torus (while other choices 
are of course possible). The mapping class group of the 
three-dimensional torus is SL{3, Z), and, as in the case of 
(2-|-l)d, is also generated by two transformations, which 
we also call the modular S and T transformations (see 
Sec. HA for details). For (3 -I- l)d topological phases 
defined on a spatial torus T^, S and T matrices can be 
introduced to describe the basis transformation of degen¬ 
erate ground states. As in (2-|-l)d, the S and T matrices 
encode the topological data of the bulk topological phase, 
such as the braiding and spin statistics of excitations. 

In (3-|-l)d, the exchange statistics of particles has to be 
either fermionic or bosonic. On the other hand, a particle 
and a loop-like excitation, or two loop-like excitations in 
the presence of an additional background loop, can have 
non-trivial braiding and can obey non-trivial statistics. 
For the Abelian topological phase described by the BF 
theory, the S matrix describes the braiding phase be¬ 
tween particle and loop excitations, while the T matrix 


has the physical meaning of a (3 -I- l)d analogue of topo¬ 
logical spins.It has been also proposed that there exist 
(3-|-l)d topological phases that are characterized by their 
three-loop braiding statistics. 

We will demonstrate that these results, obtained and 
discussed previously from the bulk point of view, can 
be obtained solely from gapless boundary held theo¬ 
ries. More specihcally, taking various examples of (3-|-l)d 
topologically ordered phases and their surface states, 
which we put on the (2 -|- l)d spacetime torus T^, we 
compute the modular S and T matrices explicitly, and 
show that they agree with the S and P matrices obtained 
from the bulk considerations. We thereby establish the 
bulk-boundary correspondence in these (3-|-l)d topolog¬ 
ically ordered phases. Along the course, we also propose 
a bulk continuum held theory which realizes non-trivial 
three-loop braiding statistics. 

N.B. Our strategy adopted in this paper is to utilize 
boundary held theories to learn about bulk excitations 
in (3-|-l)d topological phases, by establishing a bulk¬ 
boundary correspondence. One should however bear in 
mind that boundaries may have more “life” than their 
corresponding bulk, in that a given bulk topological 
phase can be consistently terminated by more than one 
boundary theory. Therefore, it would be more appropri¬ 
ate to consider a “stable equivalent class” of boundary 
theories for a given bulk theory. (See, for example. Ref. 
41.) Nevertheless, one can expect universal topological 
properties of the bulk theories may be extracted from any 
boundary theory which consistently terminates the bulk. 


A. Outline of the paper 

The rest of the paper is organized as follows. 

In Sec. H, we consider the compactified free boson the¬ 
ory in (2-1-1 )d defined on the 3d flat torus is computed. 
This (2-|-l)d theory is not necessarily tied to a particular 
(3-|-l)d bulk topological order but serves as a warm up 
for later sections. We will show its partition function is 
invariant under S'L(3,Z). 

In Sec. HI, the surface theory of the (3-|-l)d BF the¬ 
ory at level K is studied. This theory can be subjected to 
twisted boundary conditions, which are induced by intro¬ 
ducing quasi-particles in the (3-|-l)d bulk. We will show 
that the partition functions with different boundary con¬ 
ditions are transformed into each other under S'L(3,Z), 
and form a representation of S'L(3,Z). The extracted 
S and T matrices agree with the known result.We 
will also compute the thermal entropy in Sec. HID, and 
show that there is a constant negative contribution to 
the entropy. This contribution to the boundary thermal 
entropy is expected to capture the topological entangle¬ 
ment entropy defined in the corresponding (3-|-l)d bulk. 

In Sec. IV, we introduce an additional term, the ax- 
ion term or the theta term, to the (3-|-I)d BF theory. 
The theta angle has a texture (spatial inhomogeneity) 
and affects the boundary theory by twisting the quantum 
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numbers. Being static, the texture in the theta angle is 
interpreted as a topological defect, and we will show that 
the introduction of the defect makes the surface theory 
non-modular invariant, in the sense that the action of 
modular transformations is not closed within the space 
of the partition functions. 

This BF theory with the theta term motivates us to 
consider yet another theory in Sec. V, which can be con¬ 
structed by coupling two copies of the BF theory. Com¬ 
pared with the the defect system (the BF theory with 
the theta term) discussed in Sec. IV, in the coupled sys¬ 
tem, each copy can be interpreted as playing a role of a 
defect to the other. In this system, however, there is no 
externally imposed texture. We propose this continuum 
bulk field theory realizes three-loop braiding statistics 
discussed previously.^®’^^^^® On the surface, we consider 
two copies of the BF surface theories, which are coupled 
together in their zero mode sectors. We will show that, 
by computing the modular S and T matrices explicitly, 
this system exhibits three-loop braiding statistics. 

Finally, we conclude in Sec. VI. 


II. THE COMPACTIFIED FREE BOSON IN 

(2 + 1)D 

The compactified real scalar theory in (2 -|- l)d is de¬ 
scribed by the Lagrangian density 

^ . ( 3 ) 

where, for now, the spacetime is the “canonical” flat torus 
parameterized by (t, x, y). (We will consider, momen¬ 
tarily in Sec. IIB, a generic torus parameterized by six 
modular parameters.) The boson field obeys the com- 
pactification condition on a circle of radius r, i.e., 

(() = (j> + 2 tti. (4) 

This model can be exactly solved and is dual to the com¬ 
pact C/(l) gauge theory. Under the duality, the boson 
field (j) is related to the t/(l) gauge field by 

Furthermore, quantized vortices on the boson side are 
dual to quantized charges in the C/(l) gauge theory. For 
the compact U{1) gauge field theory, the monopole (in- 
stanton) proliferation leads to a confining phase and this 
process on the scalar boson side corresponds to adding a 
cos{(f)) term.^^ds process breaks the C/(l) symmetry 
in the compact boson theory, and the particle number is 
not conserved anymore. If we prohibit the monopoles, on 
the other hand, the Abelian U{1) gauge theory is stably 
gapless. 

The free boson theory can be canonically quantized: 
The corresponding Hamiltonian is 

271-fli 2 ttR 2 

H = j dx j dy [tf -b (5,,(())2 -b (dycl))'^] , (6) 

0 0 


where x and y are periodic with radius 2 nRi and 27ri?2, 
respectively, and the canonical momentum is (r := {x, y)) 

n(r) := 9t<^(r). (7) 


The canonical commutation relation is given by 




2 R 1 R 2 


E 


„*k-(r-r') 


( 8 ) 


where 5^^^(r — r') is the periodic delta function and k = 
{si/Ri, 82 /Ri) is the 2d momentum (si G Z). 

To specify the Hilbert space, we develop the mode ex¬ 
pansion of the bosonic field (j). Due to the compactifi- 
cation condition (4), the bosonic field has the following 
expansion: 


,., , rA^i , riV2 (j)o + Trot 

^ ^ Ri R 2 27r/R]i?2 


+ (poscit.r), (9) 


where N 12 € Z characterize the winding zero modes in 
the X and y direction, respectively. The Fourier decom¬ 
position of the oscillator part (poscid, r) is given by 


(poBc{t, r) 


'JR 1 R 2 ^ 2^a;(k) 

X [a(k)e-*‘^‘-‘‘''^-ba^(k)e*‘^‘+*''''’] . (10) 


According to Eq. (8), a(k) satishes the canonical commu¬ 
tation relation 


[a(k),a'f(k')] = i5k,k', (H) 

where a;(k) is the dispersion of the free boson on a Eu¬ 
clidean three-torus and given by 



( 12 ) 


On the other hand, the zero mode part satisfies 

[(('o,7ro] = 27r^f. (13) 


Owing to the 27rr periodicity of 0(t, r), the eigenvalues of 
TTo needs to be quantized according to 


TTO 


ttVq 

r-\/i?ii?2 


Mq g z. 


(14) 


To summarize, the boson field r) can be mode- 
expanded as 


Pit, r) 


(po Npt ^ riV2 ^^ 

2Try^Rjh 2ri?ii?2 i?E ^ 

-I- ^ ^ 

y'RiR2 2^/i4k) 

X [a(k)e-*“’*-*'''^ -b a1'(k)e“*+*‘''^] . (15) 
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The Hilbert space 'Ho consists of, for each winding sector 
specified by Ni^ 2 , the zero mode part and the bosonic 
Fock space for each k ^ 0. States in the zero mode 
part are labeled by the eigenvalues of tto, and hence by 
A^o- Furthermore, different winding sectors are summed 
over. In the following, the part of the partition function 
associated to the summation over A^o,i ,2 is called the zero 
mode sector. 

A. Modular transformations on 

We now consider the theory put on a generic flat torus. 
A flat three-torus is parameterized by six real parameters, 
Ro. 1,2 and a,l3,'y. For a flat three-torus T^, the dreibein 
is given by^® 

/ Ro 0 OWl 0 0\ 

= 0 Hi 0 -a 1 0 (16) 

\ 0 0 H2 / V -7 -/3 1 / 

where Hq, Hi, and H 2 are the radii for the directions r, 
x, and y, and a, /3, and 7 describe the angles between 
directions r and x, x and y, and r and y, respectively. 
The Euclidean metric is then given by 

A Be 

3fiy — 6 v^AB 

Hq -I- a^R\ + 7 ^H| —aR\ + /37H| — 7 H 2 \ 
-aRl + p-fRl Rl + (3'^Rl -jSRl , 
-7Hi -/?Hi Rl j 

(17) 

The group SL{3,Z) is generated by two transforma¬ 
tions: 

/ 0 0 1 \ / 1 1 0 \ 
t/i = 1 0 0 , 1 / 2 = 010 . (18) 

\oio/ yooiy 

Under the U 2 transformation, the metric is trans¬ 
formed as 

9^.^ {U 29 UJ)^, (19) 

which corresponds to the changes 

a—>-0 —1, 7—)-7-|-/3, (20) 

while Hq, Hi, H 2 , and /? are unchanged. 


On the other hand, Ui can be decomposed as 

/ 0 - 10 \ / 100 \ 

Ui = U[M, t/( = 1 0 0 M = 0 0 -1 

yooiy yoioy 

( 21 ) 

where U[ corresponds to the 90° rotation in the t — x 
plane and M is the 90° rotation in the x — y plane. The 
generator U[ acts on the metric as 

9^.u^ {U[9U[^)^^ (22) 

which corresponds to the changes 

Ho^-Ho/|r|, Hi^Hi|r|, n ^ 

7 —>■ —/3, ,0 —t 7 (while H 2 is unchanged), (23) 

where we have introduced 

T = a + irm, roi = Hq/Hi. (24) 

Observe also that under Hq —>■ Ho/|r| and Hi —>■ Hi|t|, 
T 2 —t T 2 /\t\^. Hence, U[ induces t —>• —l/r. 

The two transformations C/( and U 2 correspond re¬ 
spectively to modular S and T~^ transformations in the 
T—x plane, generating the SL{2, Z) subgroup of SL{‘i, Z) 
group. Combining with M, they generate the whole 
SL{3, Z) group. In the following, we call U[M as S trans¬ 
formation and U 2 as transformation. 

Moreover, the two generators of S'L(3,Z), the Ui and 
U 2 transformations defined in Eq. (18), satisfy'"*®0^ 

UiUl = Uf = H® = {UiRf = (HUi)'* = 1, (25) 

([/li?2)4 ^ ^ ^ ^ ( 26 ) 

{UiR^UifR^ = R^{UiR^Uif mod 3, (27) 

where 

R={U2 UifU^^UlU^^UiU2Ui. (28) 


B. The partition function on 

In this section, we calculate the partition function of 
the compactified free boson theory on the three-torus in 
the presence of the generic flat metric. The Euclidean 
action is given by 
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S = 


1 


271- 


(27r)2 

1 




2tvRo 2-kRi 27 rfi 2 


dr 


dx 


dy 


idr<l^r + 


2 / a^E? 


m 


-^2—I-1J {dx4'Y + ^ ^ 


+ 2 a^ {dr4>) {da:4>) + + 7 ) (5r(/<) + 2 ( ^^^^a{a/3 + 7 ) + ^/3 ) {da^cf)) {dy4>) 

Ho Hq \ 


/ R 1 R 2 

\~W 


Ih 

Ri 


(29) 


where 0 < < 2 tt are angular variables and we noted 

\/\g\ = RqRiR 2 , t = RqO^, X = Ri 6 ^ and y = 

The inverse metric 5^^ {e*^) is the in¬ 

verse of e^) is given by 


/ 


g>^- = 


1 

Ct 


a/S+'Y 

Rl 

Rl 


Rl 

a 


1 

a{af3+'Y) 1 

R'o 

Rl + 

Rl 

d2 I 

afS+o 

a{af)+^) 


(a/3+7)^ , /3= 

Pi 

HI 

+ w. 

~r p2 

Hq 


\ 




1 

Ml 


(30) 

In the operator formalism, the partition function cor¬ 
responding to the action (29) is given by the trace of the 
thermal density matrix exp(—27ri?o-ff0 over the Hilbert 
space Ro'- 


(where i,j are not summed in RiRjg'^^). I.e., 


H' = 


dxdy 
aRi 


W 


{dx(t) + P^dy(l)f + {dy(l)Y 

Hi 


R 2 


+ 2i{al3 + 7 ) 

Hq Ho 


(35) 

iQ rxo 

The mode expansion for the bosonic field (j) is still given 
by Eq. (15), where the energy spectrum a;(k) is now given 

by 


Z = Tr 


'Hq 


^ — 2t:RqH' 


(31) 


The “boosted” Hamiltonian H' and the (untwisted) 
Hilbert space Ro are specified as follows. The boosted 
Hamiltonian H' consists of the “unboosted” Hamiltonian 
H (with a = 7 = 0), and the momentum Px,y, which in¬ 
duces the boost in x and y directions, respectively: 


H' = H + i^Px + i{al3 + -r)^Py 


Ro 


Ro 


(32) 


a;(k) = ij G^^kikj 



= \ 9 ^^s^Sj = 


■ (36) 


where 


^7 = y i = x,y. 


and is defined as 

ffl2 


g^l = 


921 922 





The Hilbert space Ro is given as a direct product of the 
(33) bosonic Fock spaces each built out of a given zero mode 
state specified by iVop^. 


Next, we proceed to compute the partition function 
and study its properties under modular transformations 
of the three-torus. The partition function can be split 
into the zero mode part, which we call Zq, and the oscilla- 
(34) tor part, which we call Zqsc- Tbe total partition function 
is Z = ZoZqqq. 
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The partition function of the zero mode part is 


Zo= ^ exp 

-^0 , 1 , 2 £ ^ 


^^2 ^,2 , « A ^^2 ^TTT RoRi^^ 


- 27rr^i?2T2(iVi + ^^¥ 2 )^ - 


2 r 2 i ?2 

2TriTiNo{Ni + PN2) + 2TTijNoN2 


i?2 


(37) 


where we recall T 2 = Rq/Ri, ti = a and r = ti + zt 2 . 

On the other hand, for the oscillator part, the Hamiltonian is 


Ksc = E 
MO 


. Ri 


R 2 


w(k) + ia-—ki + i{aP + l) — k 2 
Hq ito 


a'^(k)a(k) + Eq, 


(38) 


where Eq is the ground state energy and needs to be properly regularized: 

1 


^^0= E 


sez 2 /(o.o) 


sez 2 /(o,o) 


1 


(39) 


The partition function of the oscillator part can be decomposed into the product of the partition functions of one¬ 
dimensional non-compact bosons with “mass” given by S 2 . When the “mass” S 2 = 0, 


1 _ g-2iri?o(<^(k)-|-ia^) 

-1 

1 

L 


V(t) 


= e- 2 ^-RoSo(s 2 =o) 

Si 

where rj(T) is the Dedekind eta function and r is the 2-dimensional modular parameter: 

00 

77(r):=e^ 9 := 

71 — 1 

On the other hand, the other massive part equals to 

^ g-27r_Ro-Eo(s2#0) JJ" ^ _ ^-2-KRo(u>(W)+ia^+i(aP+'^)^) 

S2/0 ,SiGZ 

= n 

S2&+ ^ ^ 

where the massive theta function 0 [,Ss 2 , 7 S 21 (d ;ij® 2 ) is defined as 

n ^ ^—2-kt 2 yj m^ + (n+a)^+27rzri {n-\-a)-\-27rib 


(40) 


(41) 


(42) 


(43) 


where 


1 _ 1 r°° 

A(m, a) = - E + (» + «)^ ~ / dk{m‘^ + k'^y^'^ 

^ ^ J — 00 


Thus the partition function for the oscillator part equals to 

2 

I 

Z‘ OSC ^S2—0^32^0 


1 

tj{t) 


n e-' 

S 2 ez+ 


Ri 

[/3S2.7S2] I 7 ^ 2 ^^ 


Together with (37), we have completed the calculation of the total partition function, Z = ZqZosc- 


(44) 


(45) 


It is instructive to compare the above partition func- ified free boson. Performing dimensional reduction, by 
tion with the partition function of the (l-|-l)d compact- 
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taking i ?2 = 1, N 2 = 0 and S 2 = 0, the partition function 
reduces to 


Z = 


1 

|?7(t)|2 


^ exp( 

7Vo,iGZ 
f 2TriTiNoNi ). 


P^o-27rrV2iVf 


(46) 


It is also straightforward to show that the oscillator part 
is invariant under M transformation and thus the total 
partition function is invariant under M transformation. 

Finally, it is also easy to check that the partition func¬ 
tion is invariant under U 2 transformation. Hence the 
partition function is invariant under the SL{i,'L) trans¬ 
formation. 


This is the partition function for the compactihed free 

boson in (1 -f l)d. III. THE SURFACE THEORY OF THE (3-|-l)D 

BE THEORY 


1. Modular invariance 


A. Bulk and surface theories 


We now show that the total partition function is in¬ 
variant under the 5'L(3,Z) transformations. 

For Zq, under I/( transformation, by using the Poisson 
resummation formula twice, we have 

Zo{t) Zo(-l/r) = \t\Zo{t) (47) 

where the Poisson resummation formula is 

(48) 

For Zosc part, under U[ transformation, the massless 
S 2 = 0 component will contribute a 1/|t| prefactor. The 
massive part is invariant under transformation, since 
the massive theta function satisfies 


a. The bulk field theory The (3+l)-dimensional one 
component BF theory is described by the action 


^bulk — 


/M 


K 

— & A da — a A Jgp — b /\ Jqy 

ZTT 




(52) 


where /r, ... = 0,..., 3, a = a^dx^ and b = 
{l/2)bp,vdx^dx'' are one and two form gauge helds; M. 
is the bulk spacetime manifold. The “level” K is an inte¬ 
ger. The three form Jgp and two form Jgy represent cur¬ 
rents of zero-dimensional (point-like) quasi-particles and 
one-dimensional quasi-vortex lines, respectively. The BF 
theory furnishes the following (bulk) equations of motion 


= 0 [b.-a] (49) 

Thus the total partition function is invariant under U[ 
transformation. 

Under the M transformation which is basically a 7r/2 
rotation in the x — y plane, the partition function for the 
zero mode part becomes 


Zo A Zo = 




exp 


Nn 


T^Rq 

2 r 2 i?ii ?2 




—2TTr^RoRiR2 





— 2TriaRo[l3Ni — N2]No — 2nijNiNo^ ■ 

(50) 


do — jgvi r, db — Jgp- (^^) 

Z7T ZTT 

The BF theory implements a non-trivial fractional 
statistics between quasiparticles and quasivortices. To 
see this, we consider the following conhguration of quasi¬ 
particles and quasivortices: 

Jgp = 5^{C), Jgy=52{S). (54) 

Here C and S represent the one-dimensional wold-line and 
the two-dimensional world-sheet of quasiparticles and 
quasivortices, respectively; do-ni-N') is the delta func¬ 
tion {D — n)-form associated a submanifold Af C M, 
where D — n = dimAI — dim A/”. By definition, for any 
n-form 




di)—n(A7) A 


Hence, for example. 


(55) 


Therefore the invariance of the zero mode part of the 
partition function can be seen from relabeling. 



Jqp O — 




A6 = 



(56) 


[iVi r-iV 2 

M ■. In2 ^ < iVi 
[No [ao 


Some useful properties of the delta function forms are 
summarized in Appendix B. 

In the presence of these quasiparticles and quasivor¬ 
tices, we now integrate over a and b to derive the effective 










action for Jqp and Jg„. Since the theory is quadratic, this 
can be done by solving the equations of motion. These 
equations, up to a closed form, are solved by 

Stt Stt 

h — d Jqp^ d Jqv 

(If the spacetime is trivial, by the Poincare lemma, a 
closed form is exact. If so, such exact term does not affect 
our final result since, for an arbitrary closed submanifold 
Af, f S{N'){d4‘) f = 0.) From the formula 

(B4), d~^Jqp and d~^Jqy are determined as 


b. The boundary theory On a closed manifold, the 
BF theory is invariant under gauge transformations a —> 
a + dip, where ip is zero form, and 6 —>■ b + d(^, where ^ is 
one form. In the presence of a boundary (surface), there 
may appear gapless degrees of freedom localized on the 
surface. The action describing the boundary degrees of 
freedom can be inferred by adopting the temporal gauge 
flo = biQ = boi = 0 (i = 1,2,3), solving the Gauss law 
constraints e^'^^^dkbij = = 0 by = dkP, 

bij = diQ — djC,i, and then plugging these back to the 
action. The resulting action is^®^"'® 


d ^Jqp = 62(T>), where &D=C, 

d~^ Jqv = di(y), where = 5. (58) 


where the two-dimensional manifold T) and the three- 
dimensional manifold V are chosen such that dT> = C 
and dV = S. They are not unique, but different choices 
lead to different d~^Jqp^qv which differ by closed forms. 

Substituting these solution into the action, 


27r f 

Sbulk — ^ / (^ Jqv) Jqp 


= -^LkiS,C), 


(59) 


where Lk is the linking number between Jqp and Jqy. 
Hence, 


J V[a,b]e^^>’-‘'^ ^^-^Lk(s,c) 

In the last line, we assume the world-line C consists of tra¬ 
jectories Ci of many quasiparticles each carrying charge 
Qi e Z: Jqp = 63 ( 0 ) = Y.idih{Ci)- Similarly, the world¬ 
line S consists of trajectories Si of many quasivortices 
each carrying charge \i G Z: Jqy = 52{S) = XiS 2 {Si). 
The fractional phase (when |K| > 1) in Eq. (60) repre¬ 
sents statistical interactions between quasiparticles and 
quasivortices. 

Once the coupling of the gauge fields to the currents 
is prescribed, it also specifies the set of Wilson loops and 
Wilson surfaces included in theory (see Eq. (56)). If the 
theory is canonically quantized on A4 = R x S, the set 
of the Wilson loop and Wilson surface operators of our 
interest is 


exp im 



expin 



(61) 


where m, n are integers, and L and S are arbitrary closed 
loops and surfaces in S. These operators satisfy the com¬ 
mutation relations. 


2Tri 




^imf^ a^inf^ I 


(62) 


where I{L,S) is the intersection number of the loop L 
and the surface S. 


SdM 



dtdxdy 


K 

27r 


- V{p,0 


(63) 


where i,j = 1,2. Here we have added the potential 
V(p, ()), which originates from microscopic details of the 
boundary and is non-universal. This boundary action 
can be obtained from the the free scalar and the t/(l) 
Maxwell theories by imposing a self-dual (or an anti-self- 
dual) constraint, e>^''^dvC,x = 

The appearance of the gapless degrees of freedom on 
the surface deserves more comments. In particular they 
should be contrasted with the gapless edge theory of the 
(2-|-l)-dimensional Chern-Simons theory. For the single¬ 
component Chern-Simons theory in (2 -|- 1) dimensions, 
the boundary is described by the single-component chiral 
boson theory, which is stable and cannot be gapped out. 
The appearance of the gapless edge theory is necessary 
since the bulk theory is anomalous and the anomaly must 
be compensated by the degrees of freedom living on the 
edge. 

On the other hand, the surface theory of the BF the¬ 
ory in (3-1-1) dimensions (also in (2-1-1) dimensions) can 
be gapped out by adding suitable perturbations (if we 
do not require any symmetry). In other words, there is 
no anomaly protecting the gapless nature of the surface 
theory. Nevertheless, the appearance of the surface the¬ 
ory (63) can still be understood in terms of an anomaly. 
While the BF theory (both in (2-1-1) and (3-1-1) dimen¬ 
sions) is equivalent to the topological Zk gauge theory, 
which does not have gauge anomaly on a manifold with 
boundary, the continuum action of the BF theory ar¬ 
tificially preserves t/(l) symmetry. The BF theory is 
anomalous under C/(l) in the presence of a boundary, 
and this anomaly must be canceled by gapless degrees 
of freedom living on the boundary. As in the bulk the¬ 
ory, the boundary theory is invariant under the artificial 
U{1) symmetry; this symmetry translates the boson field 
p ^ p + const. If the 17(1) symmetry is strictly pre¬ 
served, the gapless boundary theory cannot be gapped, 
as can be inferred easily from the fact that any gapping 
term of cosine type cos(np + a) violates the 17(1) symme¬ 
try. (If we use the dual picture of the compactified boson, 
i.e., the compact 17(1) gauge theory, the 17(1) symmetry 
is equivalent to prohibiting monopoles.) On the other 
hand, once we relax the 17(1) symmetry, which is, from 
our point of view, an artificial symmetry after all, this 
gapless boundary theory can be easily gapped out by 
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adding some relevant perturbation and is not stable at 
all. 

While this gapless surface theory is not stable at all, 
it does encode topological data of the bulk, as we will 
demonstrate later. Let us for now discuss, in more de¬ 
tail, the connection between the bulk excitations and the 
fields living on the boundary. In the following we choose 
= 5"^ X E, where the spatial manifold E is a solid 
torus, E = X , and hence dM. = . Let us first 

consider a quasiparticle current consisting of a quasipar¬ 
ticle carrying no units of charges (no G Z): 

jqp{x) = ^0 ^ - X{t)] (64) 

where L is the world-line of the quasiparticle, and the co¬ 
ordinate X^{t) represents the trajectory of the particle. 
For the quasiparticle at rest, X^’^’^(r) = const. = 

j°p{x) = noS^^\x-X). (65) 

Integrating the equation of motion over the total space, 

^ J d^xe°’-^''dibjk = J d^xjgp = no. (66) 


Using Stokes’ theorem, db = fg^ b = 
(1/2) Jg^bijc'^^d'^x, and substituting 6^ = diQ — djQ, 
this reduces to 

[ d^xeijdiCj =‘^no. (67) 

JO'S ^ 

Hence adding a quasiparticle in the bulk corresponds to 
introducing flux on the surface. 

Similarly, let us consider to introduce a quasivortex 
source: 



0/3 


dX^{a)dX^{a) ^f,, 
da°‘ dal^ ^ 


X{a)], 

( 68 ) 


where S is the world-surface of the quasivortex, and the 
coordinate X^{a) represents the trajectory of the par¬ 
ticle in spacetime, and n is an integer. Let us consider 
a straight quasivortex at rest, stretching along a non- 
contractible cycle of the bulk solid torus. For conve¬ 
nience, this direction is taken as the a:-direction (Fig. 
1). Then, = 0 and 

jll{x)=n^5{x^-X^)5{x^-X^). (69) 

where X^’^(cr) = const, and we have renamed the integer 
n as 712. Integrating the equation of motion over space, 

Li X ^ [ dydz diUj = n2 x Li, (70) 
Stt j 

where 7,/ = 1,2 and Li = 27ri?i is the length of the 
quasivortex stretching in the x-direction, and we noted 



FIG. 1. (a) The presence of a point-like quasiparticle in 

the bulk (solid torus) induces a fractional flux on the spatial 
boundary S (torus), (b) The presence of a quasivortex line 
in the bulk twists the boundary conditions of the surface the¬ 
ory. Here and in (c), the bulk is presented as a filled cylinder 
where the top and the bottom of the cylinder are identified. 
The shaded surface is a sheet of the branch cut which em¬ 
anates from the quasivortex, and intersects with the spatial 
boundary (depicted by a wavy line). The surface excitations 
experience a twisted boundary condition as they go through 
the branch cut. (c) Similar to (b), a bulk quasivortex, which 
creates a branch cut on the surface which now goes along a 
different cycle of the surface, is depicted. 


the flux Oj is independent of xx. Using Stokes’ 

theorem, and substituting Oi = di^p, 


j) dydyp = —712. (71) 

Hence introducing a quasivortex (quasivortices) along the 
non-contractible loop in the bulk corresponds to intro¬ 
ducing winding of the scalar boson on the surface. 

One may wish to develop a similar argument for a 
quasivortex (quasivortices) stretching in the j/-direction. 
(Fig. 1 (c)). It should however be noted that once we 
fix our geometry as above (Fig. 1 (b)), loops running 
in the y-direction are contractible in the bulk. In other 
words, if one constructs a solid torus by filling “inside” of 
a two-dimensional torus, one needs to specify one of non- 
contractible cycles on the two-dimensional torus, such 
that after filling, this cycle now is contractible in the 
sold torus. 


B. The surface theory and quantization 

We now proceed to the canonical quantization of the 
surface theory. We start from the surface Lagrangian 
density 


C = 


K 

27r 


{€,jdiCj){dt(p) 





^2 


1 

2 A 2 


G^^dipdjp. 


(72) 


The boson field (p is compact and satisfy 


p = p + 2 'k. 


(73) 
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I.e., physical observables are made of bosonic exponents 

exp[iTO(/?(t, r)], rnGZ, (74) 

and the derivative of the boson fields (current operators). 
The winding number of ip is quantized, in the absence of 
bulk quasiparticles, according to 

dx^diip = 2TTNi, Ni £ Z, (75) 

where i = 1,2 and i is not summed on the right hand 
side. On the other hand, the gauge field Q is compact, 
meaning that physical observables are Wilson loops, 

expim / m £ Z, (76) 

Jc 

where C is a closed loop on 9E = . The flux associated 

to Q is quantized, in the absence of bulk quasiparticles, 
according to 

dxdy CijdiQ = 2 ttNo (77) 

where Nq is an integer. The canonical commutation re¬ 
lation is 

o jr7 

[ip{t,r),e^jdiQ{t,r')] = (78) 

In the following, we fix Ai and A 2 according to 


where the eigenvalues of /3o and /?i ,2 describes the flux 
(associated with the gauge field a in the bulk), and the 
winding of the ip field, respectively. The quantization 
conditions of these variables will be discussed momen¬ 
tarily. Reflecting the compact nature of the p and C,j 
fields, the zero modes are compact variable -|-27r 

(/i = 0,1,2); For ag, the compactification condition 
comes from the fact that physical observables are given 
as bosonic exponents (74). Similarly, for 01 , 2 , that phys¬ 
ical observables are given in terms of Wilson loops (76), 
and that these Wilson loop operators must be invariant 
under large gauge transformations imposes the compact¬ 
ification condition, 01.2 = 01,2 + Stt. 

From the commutator [p{t, r), eijdiQ{t, r')], we read off 

[a(k),a'''(k')] = i5k,k', 

= (82) 

From the compactification condition 00 = 00 + 27r, /?o is 
quantized according to 

/3o = ^, MoSZ. (83) 

This quantization condition translates into 

J dxdy Cij d, Q = . (84) 


( 2 ^)^ 

K 2 A 1 A 2 


(79) 


This choice is convenient since it gives rise to the same en¬ 
ergy dispersion as the compactified free boson discussed 
in the previous section. 

To proceed, we consider the mode expansion of the 
fields. The equations of motion are 


-^eijdACj + ^G"^d,djp = 0 , 

ZTT A2 

-^eikdidtp + -^eikdiieijdiCj) = 0. (80) 

zTT Ai 


The mode expansion consistent with the equations of mo¬ 
tion are 


+(r) = «o + ^ + ^ 

ill 1 x 2 


1 


1 


■E: 


x/HTH; V 2K2A, i..(k)i/2 


X [a(k)- 


— ik-i 


t(k).+*k-rl 


0(0 = 


—^ H-^ 

27ri?, 27ri?ii?2 


-x5. 


4.2 


1 


4 E 


x/Rjh V 87 r 2 a;(k) 3/2 1 


X [a(k). 


— ik-i 


t(k).+*k-rl 


(81) 


Compared with the quantization condition (77), the flux 
is now quantized in the fractional unit. We will separate 
Mq into its non-fractional and fractional parts as 

Mq = KiVo + no, Nq £Z, no = 0,..., K — 1, (85) 
and write the quantization condition of /3o as 

/3q = Nq + uq/K. ( 86 ) 

The quantization condition of /3i_2 can be discussed 
similarly. From the commutator [p{t, r),eijdiQ(t, r')], we 
infer 

O'JT'j 

hjdMt,r),Q{t,r')] = (87) 

which implies 

077*7 1 

[I3ua2]-[l32,ax]=- — -. ( 88 ) 

One can choose, for example, 

O'TT'} T 

[/3i,a2]=0, [/32,ai] = ^-. (89) 

This choice may be consistent with the previous consid¬ 
eration from the bulk point of view, and in particular 
with the comment below (71). I.e., this choice may cor¬ 
respond to choosing which non-contractible loops on the 
surface are contractible in the bulk, when forming a solid 
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torus starting from the two-dimensional torus by filling 
its “inside”. 

From the compactness of the gauge field Q, the zero 
modes satisfy ai = ai + 2 'k, which imposes the quantiza¬ 
tion condition 

M 9 

P 2 = M 2 G Z. (90) 

K 

As before, we split M 2 into the fractional and non- 
fractional parts, 

M 2 =KN 2 +n 2 , N2e1, 712 =0,...,K-1. (91) 

With this, the boson field obeys the twisted boundary 
condition 

(p{t, x,y + 27ri?2) = (fit, X, y) + 2 tt (^N 2 + . (92) 

While above consideration allows winding in the y- 
direction but not in the x-direction, in computing the 
partition functions of the surface theory in the next sec¬ 
tion, we consider winding in both directions, 

(p{t, X + 2TTRi,y) = X, y) + 2 tt (^Ni -h , 

(p{t, x,y + 27ri?2) = (f{t, X, y) + 2 tt (^N 2 + , (93) 

That is 


where $ is a single component boson theory compactified 
as $ = $ -I- 27r, and obeys the canonical commutation 
relation [$(a;), i9a;<I)(x')] = (2Tri/K)S(x — x'). The zero 
mode part of <&, defined by the mode expansion 

Ht, x) = ^o-pit + x)+iJ2 , (97) 

n/0 

satisfies [d>o,p] = i/K- This then suggests the quantiza¬ 
tion rule, p = (integer)/K, and the boundary condition 
of the chiral boson field 

#((,* + 2 ,) = <!.((,08) 

Thus, the canonical quantization naturally leads to the 
twisted boundary condition of the chiral boson field. 

Quantization of the surface theory with the above 
twisted boundary conditions gives the spectrum of lo¬ 
cal as well as nonlocal (quasiparticles) excitations, which 
obey untwisted and twisted boundary conditions, respec¬ 
tively. Once we specify the boundary condition (with 
some integer vector 77p=o,i,2)5 the theory is quantized 
within one sector (labeled by the equivalence class [n] 
with the relation n = n + KA where A is a vector with 
integer entries) of the original spectrum. For this surface 
theory, there are sectors in this compactified theory 
and is consistent with the ground states of single com¬ 
ponent BF theory defined on T^. 


A=1,2 — A^1,2 + ^1,2/K- 


(94) 


To summarize, in the presence of twisted boundary 
conditions, the mode expansion of the fields are given by 


diPit, r) = 


N,. 


h/K 


Ri 


—i 


'JR 1 R 2 V 2K2Ai ^w(k) 


E 


ki 


X _a(k)e-“('')‘-*‘''^ -k a^(k)e+*‘^(‘')‘+*''''' 


C. The partition functions and modular 
transformations 

Now we compute the partition function (coupled to the 
metric): 


= Tr-w 

n.nQn-^n2 




(99) 


where Rnonin -2 is the Hilbert space twisted by no,ni,n 2 
fractional quantum numbers, and 




Nq + TZq/K 
27ri?ii?2 


—i 

\/ R 1 R 2 



E 

k/O 


X 


a(k)e-*“'(k)‘-*k-r 


„t(k)e+“(k)‘+*k-r 


(95) 


The above consideration is somewhat analogous to the 
quantization of the chiral boson theory that appears at 
the edge of the (2-|-l)d Chern-Simons theory at level K. 
The (l-l-l)d chiral boson theory defined on a spatial circle 
of radius 27r is described by the Lagrangian density 


H' = H 


.aRi 
■ 1 

E? ^ 
Kq 

f K^A 
^ = / 'i^dy 


*(a/3 + 7)^ 

Ho 

477^ 

J^2;^2 Hb'C'iQ , 


f K 

Pi= dxdy—{€imdiCm){dtip). 


Pyi 

dipdjip , 


( 100 ) 


The calculation of the partition function goes in paral¬ 
lel with the calculation presented in the previous section 
for the free boson theory. To see this, we note, from the 
equation of motion. 


GjdiCj 


KAi 

27r 


dtP 


( 101 ) 


C = 


K 

47r 






(96) 


up to a constant term. Thus, in terms of the Hamil¬ 
tonian density and the commutation relation are given 
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by 

Kt,r),9*<^(t,r0] = ^<5(^)(r-r0. (102) 

By introducing the rescaled field, 

^ (103) 


the Hamiltonian and the commutation relation can be 
made isomorphic to those of the free boson theory. The 
compactification condition of the rescaled boson field is 

r = K^. (104) 

The partition function i^g computed 

from the partition function of the free boson theory. The 
zero mode part of the partition function for each excita¬ 
tion sector is obtained from Zq (Eq. (37)) by making re¬ 
placement TVo “t KNo+uq and Ni —>■ iVi-|-ni/K {i = 1, 2): 


^non\n2 


-^0,1,2 £ ^ 


7rK^T2 

2r2i?2 




+ 27rjTiKiVo 


Ni + 02 


where we have introduced the notation 


iV^ := + n^/K. (106) 

For the oscillator part, since the Hamiltonian is the same 
as the oscillator part for the compact boson, the parti¬ 
tion function is exactly the same as the free boson case 
presented above. Thus we have 


^S2=0'^S2#0 


1 

ry(r) 


n 


S2GZ+ 


0 


-1 

[/3s2,7S2] 



(107) 


The total partition function for each sector is = 

Znon^u^Zosc■ 

Although the surface theory of the (3-|-l)d BE theory 
resembles the compactified free boson discussed in the 
previous section, these theories are physically different. 


2r^7ri?2'r2 


A^i + /3A^2 


2 2r‘^TTRoRi ~,2 


i?2 


-m 


2Tri 


*7K7Vo7V2}, 


(105) 


For the compactified boson, the partition function is in¬ 
variant under the S and T modular transformations: It 
is anomaly-free and a well-defined theory on the (2 -|- l)d 
spacetime torus. On the other hand, for the surface the¬ 
ory, the partition function for each sector is not modular 
invariant and thus it is not a well-defined theory on the 
(2-1-l)d torus. It should be regarded as the boundary the¬ 
ory of a higher-dimensional topological phase. There are 
sectors determined by three quantum number nQ i 2 
and they form a complete basis under S and T modular 
transformations, as we will show now. 

Under M transformation, quantum numbers are trans¬ 
formed as 


M:In2 + ^ ^ < A^i + ^ (108) 

[No + ^ [iVo + t 

To discuss U[ transformation, we use the Poisson re¬ 
summation to rewrite the summation over Nq and A^i in 
rewrite the zero-mode partition function as 


^nor 


1 


K| 


E 


exp< - 


27riTiK Ml 
|r|2 IT 


2r2i?. 


^^-2 ,.2 2 r 27 ri ? 2 r 2 

'R20 ' |t|2 


Mo 

K 


- 7^2 


27rr^i?oi?i 

i?2 




Mo . ' 

— - 7iV2 


-h 2TTi02Mi + ^^^Mi 

K 


27rmo 

K 


Mn 


(109) 


Let us introduce 


Ml := KNii + n'o^ Mq := KN[ + n'l 


( 110 ) 
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Then, the partition function can be written as 
1 


^nor 


K\t 


E E 

N2,Mq^G'^ 

27rzriK 


TTT2K‘^ ~,2 2r'^TTR2T2 

2 r 2 i? 2 |T |2 0 |r|2 


12 ^0 


N[--^N2 + 27riK/3iV2iV^ ■ 


N[ - 77 V 2 

27rzninQ 27rinoni 


2 2Trr^RoRi 


i?2 




K 


K 


( 111 ) 


where X^neZK X]n=o- From these expressions, under 
U[ transformation. 


^jjr^2^'jnonin2 _ }_ 'y g^(no"'i+"i"o)2^"o"i"2 
njj i^Zk 


( 112 ) 


Here, we made a convenient choice Ai = 1/K, i.e., 2r^ = 
K. This is the same as the character of the edge theory 
of the (2+l)d Zk gauge theory in its topological phase. 
The effective Lagrangian density of the edge CFT is by 

C = (115) 


Combined with the M transformation, we can write down 
the modular S and T matrices: 


S ,——(1 ,p-¥(">2-noO 

T >=5 >5 >6 11131 

>ni,n^ ^no ,^0 ^ni ,nj^ ^n2 ,^2 


This result is consistent with previous works. Refs. 33 and 
34, and also 36, where the action of the modular trans¬ 
formations are calculated in the bulk. (See also other 
related works: Refs. 35, 37-39.) In terms of the bulk 
physics, the S matrix describes the braiding phase be¬ 
tween particle and loop excitations, whereas the T ma¬ 
trix encodes information related to (3 -I- l)d analogue of 
topological spins.(See also Refs. 35, 37-39.) The ex¬ 
act agreement between the S and R matrices calculated 
in the bulk and the boundary suggests there is one-to- 
one correspondence, the bulk-boundary correspondence 
in (3-hl)d. 

The compnted S and T matrices (113) are expected to 
be consistent with the algebraic relations in Eq. (27): 
As in (l-l-l)d CFTs, together with the charge conju¬ 
gation matrix C, S and T matrices should obey essen¬ 
tially the same algebraic relations as Eq. (27). Assum¬ 
ing the charge conjugation matrix is unity, C = 1, we 
have checked, for the case of K = 2, 3,4, 5, the S and T 
matrices satisfy all the above constraints except the last 
equation in Eq. (27). 

Before we leave this section, as we have done in the 
previous section, it is instructive to dimensionally reduce 
the partition functions of the surface theory of the (3-|-l)d 
BF theory. For each given sector, after dimensional re¬ 
duction, the partition function is given by 


= S e>^p{-»r,K(jV„ + |)’ 

^ ' Afo.iez ^ 

- 7rT2K (iVi + 

+2^iriK (iVo + (iVi + ^) } • (114) 


where K = and V is a symmetric and positive defi¬ 
nite matrix that accounts for the interaction on the edge 
and is non-universal. The characters defined in Eq. (114) 
can be simplified as 


Xabir) = 


\viT-)\ 


■ \ ' /74i< 


2 / 

s.t 


(116) 


where a = no and b = ni. There are characters in 
total. Under the S and T modular transformations, they 
are transformed as 

Xabir + l) = e2”Trxab(T), 

Xab(-l/r) = ^ E Xa'b'iT)e-^^^'^-^. (117) 

a' 


D. Entropy of the boundary theory 

In this section, we compute the thermal entropy 

St ^ [TIn xa], (118) 

obtained from the partition functions of the boundary 
theory discussed above. Here, Xa is the partition function 
in the sector labeled by a = (no, 711 , 712 ), and 

l/r = 27rRo (119) 

is the inverse temperature. 

While St is defined for a system with a real (physi¬ 
cal) boundary, it is expected to carry information on the 
universal topological part of the entanglement entropy 
(the topological entanglement entropy). The latter is de¬ 
fined for the bulk system (the BF theory) defined on a 
manifold without a physical boundary, and obtained by 
integrating out (tracing over) a subregion B (compliment 
to, say, subregion A).®49-5i 
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We are interested in the entropy St in the limit 
Ri/Ro —t 00 and R 1 /R 2 00 . (We could also equiva¬ 

lently take the limit with i?i and i ?2 exchanged, in which 
case, we have to resum differently but the result would 
be the same.) To evaluate the entropy in this limit, we 
first make use of the 5-modular transformation,^^ and 
write 

5T = ^[Tln(5^X,(-l/r))]. (120) 

In the above limit, only the identity charac¬ 
ter gives rise to the dominant contribution, 
linifii/iJo^oo,fli/fl 2 ^ooXa(T) = 5°xo(-l/'r), as seen 
from Eq. (111). Hence 

St\ri/Ro^od,Ri/R2^oo = [^ la ( 5 ° xo )] ■ ( 121 ) 

Then using the modular S matrix computed in the pre¬ 
vious section, 

^t\ri/Rq — ^oo,Ri/R 2 — ^00 (RP 

= -lnK+^[Tlnxo]. (122) 

The first term is the subleading term, and identical 
to the bulk topological entanglement entropy, although 
St and the entanglement entropy are defined differently. 
The second term is the extensive piece, which basically 
corresponds to the entropy of the free boson and is the 
usual leading order term. (When St is interpreted as 
the entanglement entropy, the second term corresponds 
to the area law term.) 

IV. THE SURFACE THEORY OF THE (3 + 1)D 
BE THEORY WITH THE 0 TERM 

Recall that in the surface theory of the BF theory dis¬ 
cussed in the previous section, there are three quantum 
numbers Mg, 1 , 2 , which we wrote in terms of their non- 
fractional and fractional parts as 

M^ = N^ + ^, ^ = 0,1,2. (123) 

These quantum numbers in the surface theory can be 
interpreted as arising from the presence of bulk quasi¬ 
particles or quasi-vortices; Mi ,2 represents the frac¬ 
tional winding number of the tp field induced by a bulk 
quaxi-vortex, whereas Mg represents a fractionalized flux 
threading the surface induced by a bulk quasi-particle. 

In this section, we consider the following “twist” of the 
quantum number 



in the surface theory of the BF theory, where Qi ,2 are 
fixed integers, and we have introduced the notation 

Q X M ■-Q1M2-Q2M1. (125) 

This twist can be induced by considering a modification 
of the BF theory by introducing the 0 term (axion term). 
In the next section, we will consider a similar twist to 
discuss three-loop braiding statistics. 


A. The BF theory with the 0-term in (3-|-l)d 


We motivate the twist (124) by considering the fol¬ 
lowing modification of the bulk BF theory by adding a 
0 -term: 


Shulk 



— b A da -^d0 Aa Ada 

2 tt Stt^ 


— a A Jqp — b A Jqy 


(126) 


In the second term (the-0 term or axion term), p is a 
parameter, specific value of which will be discussed later, 
0 is a non-dynamical background field, and we consider 
an inhomogeneous but time-independent configuration of 
0, which will be specified later. Compared to the stan¬ 
dard form of the 0 term, Qda A da, we have done an 
integration by part and put the derivative acting on 0. 
Since the 0 field is non-dynamical, we will interpret the 
presence of the 0 term as an introduction of a static 
defect. In Ref. 53, a similar effective action has been 
proposed to describe the thermal and gravitational re¬ 
sponse of topological defects in superconducting topolog¬ 
ical insulators.We also note that the BF theory with 
the 0-term, Qda A da, has been proposed to describe 
the fermonic and bosonic topological insulators. In Ref. 
37, the BF theory with the 0 term was used to discuss 
three-loop braiding processes. 

To see the 0-term induces the twist (124), we assume 
the following configuration of the 0-field: 


e{x,y,z) 


Qix Q 2 y 
i?i i ?2 ■ 


(127) 


where Qi ,2 are fixed integers. From the equation of mo¬ 
tion. 


(128) 

By plugging the first equation into the second, these 
equations of motion reduce to 


MoMo + 


Q 1 M 2 - Q 2 M 1 


K 


= Mo + 


QxM 


(124) 


K 

27r 

K 

47r 




e^''^^d, 2 bxp = - 


P 

27rK 


K 


duQjqv + jqp- 


(129) 
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In the presence of quasiparticle and quasivortex sources, 
(65) and (69), the equations of motion integrated over 
space are 


^ j dydze^^^^diQj = Ml, 

^ J dxdz dittj = M 2 , 

^ ^ d^xe°^^^d,bjk = -^Q^M, + No, 


(130) 


where Mi^ 2 ,No € ^ and we noted 

dS J d^xj°l = X (27ri?,) = 27rQ,M, (131) 

(i is not summed over). These can be reduced to, by 
using Stokes’ theorem, 

f 2 tt 

J dydy^} = —iVa, 

f 2 tt 

/ dxd,^ip = z^Ni, 

J d'^xeijdiCj = + ‘^{Q X N), (132) 

where Mi = —N 2 and M 2 = Ni. Hence, upon choosing 
p = 1, in the presence of the defect field 0, the quantum 
number Nq in the surface theory is “twisted” as in Eq. 
(124). 

We observe that the following action 


SLik = I 


. a - \b - 

Jqp^fJ^ Jqv 2 


ZTTiy j 

(133) 


we note that the quantization rule of the zero modes are 
determined both by (a) the canonical commutation rela¬ 
tions and (b) the compactification conditions. The com- 
pactification condition is determined by declaring physi¬ 
cally observable Wilson loop operators. This in term is 
determined from the coupling of the theory to the cur¬ 
rent. Therefore, in the original theory, (a) is abnormal 
but (b) is normal. In the modified theory, (a) is normal 
but (b) is abnormal. In the next section, we demonstrate 
this by deriving the quantization conditions (132), de¬ 
rived from the bulk point of view here, in terms of the 
boundary theory of In Appendix A, we quantize 

the boundary theory of the original theory, Shuik, to de¬ 
rive the quantization rule (132). 


B. The surface theory and partition fnnctions 


1. The compactification conditions and quantization rules 

We now proceed to consider the surface theory of the 
bulk theory (133). Without sources, the surface theory is 
described by the same Lagrangian density as the surface 
of the BF theory, (72), and hence has the same canonical 
commutation relations. This immediately means that the 
oscillator part of the surface theory can be treated in 
exactly the same as before. On the other hand, reflecting 
the abnormal coupling of the gauge fields to the currents 
in the bulk action the compactification conditions 

of the boundary fields (p and C are modified, as we will 
now discuss. 

As we noted earlier, the coupling to the current can 
be written, e.g., J d'^x a. Thus, introducing a 

proper current corresponds to introducing a Wilson loop. 
If we now consider a Wilson line L that is spatial, and 
that ends at the boundary. 


shares the same equations of motion, Eq. (129), as the BF 
theory with 0 terms, (126). Hence, the boundary theory 
derived from has the same quantization rules of the 
zero modes as the boundary theory of S^uik- In the next 
section, we will consider the boundary theory derived 
from and its partition functions. 

To contrast the two theories Suuk and we note, 

in Sbuiki that the coupling to the currents are “normal” 
while the commutators are “abnormal”, in the sense that 
the commutators among fields a, b are modified due to the 
presence of the theta term. On the other hand, in 
the commutators are normal (the same as the ordinary 
BF theory) while the coupling to the current is “abnor¬ 
mal”. (Since that the commutators are the same as the 
ordinary BF theory, and the corresponding bound¬ 

ary theory can be analyzed in a complete parallel with 
the BF theory - a practical reason why we will consider 
on in the following - expect for the zero mode part.) 

In spite of these differences, these theories lead to the 
same quantization conditions (the same “lattice” of quan¬ 
tum numbers) of zero modes. To see how this is possible. 


a = dip = (p = (p{dL) 


(134) 


IdL 


where we noted dL is a point, and we have used the 
solution to the Gauss law constraint, Ga = da<p {a = 
1, 2, 3). Thus, 


expim J a = exp[im(p{dL)] (135) 

This means that ip is compactified with the radius 27r. 

Let us repeat the same exercise for the coupling to the 
quasi vortex current: 




h p 


's 


qv 2 L 

b — A d<d 

2 ttK 




(136) 


where S is the world surface of a quasivortex (quasivor¬ 
tices). In the presence of a boundary and using b = d^, 
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this is evaluated as 


Inverting this relation, 


Is 


dC- ^^d>fAde 
^ 2 ttK ^ 


Ids 




(137) 


where the boundary of the world sheet is on the surface. 
We thus have a Wilson line on the surface: 


exp tm 


Ids 


c- 


27 rK 


(fi A dQ 


(138) 


= Su (148) 


where 


uo = (1,0,0)^ 


, U2=(^^,0,l] . (149) 


We now consider the case where dS is along the x- or 
y- cycles. Recalling the mode expansion Eq. (81), and 
noting (p/27rK) (pdQ = (p/K)Q:oQi the zero modes 
enter into the integral ^ through the 

following combinations 

ai — ^QiciO; 0:2 — ^Q2 o;o- ( 139 ) 

Together with ao, the following three linear combina¬ 
tions 


Renaming the integers as mo —> Nq, mi —> —IV 2 , and 
m 2 -A Ni, Eq. (148) is nothing but the quantization rule 
(132). 

2. The partition functions 

With the twist (124), we can now write down the zero 
mode partition function. Let us recall the partition func¬ 
tion of the BE surface without the theta term, 
defined in Eq. (105). For later use, we write as 


a = 0,1,2 (140) 


are angular variables, where 



Noting the commutation relations among zero modes. 


Znon.n,^ ^ /k (Mq M,) 

= /K(KiVo + no,KiV, + n,), (150) 

where /k is defined by the summand in Eq. (105), and 
recall = ¥^N^+n^. We will call the partition function 
resulting from the twist . It is given by 


[ao,Po] — [ai,—/32] — [ 02 ,/3i] — (142) 

we consider the linear combinations 

^ = (/3o,-/32,/3i), a = 0,1,2, (143) 

where w°‘ are translation vectors reciprocal to v‘^: 

« = ^6 ■ (144) 

Explicitly, they are given by 

Wo = ^ 1 , , wi = ( 0 , 1 , 0 ), W2 = ( 0 , 0 , 1 ). 

(145) 

Then, in the “rotated” basis, the commutation relation 
takes the following canonical form: 


^qTqT= E /k (KiVo + no + ^^,KiV, + n,y 
Wo.i, 26 Z ^ 

(151) 

To proceed, we write 

Mi = KiVi -\- Tii = -|- Kti -|- Tii, 

Q, = KR, + r, = + Ks, + (152) 


where new integers Ni,Ri,Ri and Zk variables ti,Si,ri 
are introduced. In the following, we will show that the 
zero mode partition function depends on Qi only through 
Ti, and hence can be denoted as and that the 

partition function can be written as 

(153) 

il,2£^K 


(146) 

Due to the compacticity of w^Pi, takes on values 

WhP,. = ^ X mb, mb=i,2,3 G Z. (147) 


where we have introduced 

Rq ■= ng + s X n + r X t mod K, 

{ ho := Kuq -I- (r X n) 
hi := Kti + rii 

h2 := Kt2 + n2, 


( 154 ) 
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(i.e., ng = (ng + sxn + rx t)%K) and is defined 

by 


j^nor. 


E a( 


K/l„ + ^,K"W,+ ».). 




K 


(155) 


To show Eq. (153), we start by writing the partition 
function in terms of variables introduced in Eq. (152): 


^nQnin 2 

^Qi^K^Ri+Ksi+n 


= E E 

AfoXi, 2 GZti, 26 ZK 

X /k ^ — j + Kti + . (156) 

By further introducing 

oq = ug + s X n + r X t, 

Ao=No + {KR + s) X N + iixn + r X N, (157) 
and noting the equality 

KiV„ + „„ + 24M=K.4,+a„+Li^ (158) 


Then, 


ynQnin 2 

^Q-^K^R-+Ksi + 


K 


n= E E 

.^ 0 X 1 , 2 GZ il,2GZK 


K 


^ /k ^K^dg + og H-+ Kti + Ui ) . (159) 

We now fix ti _2 and consider 


^710711712,^1^2 

^Qi^K'^Ri+Ksi+r, 


= E 

No,N-i 2 & 


X /k ^K^dg + og H-^—, K^Xi + Kti + Ui J . (160) 

Note that once ti _2 are fixed, og is fixed. Converting the 
summation over Xg to a summation over Aq, 


^71o71i7l2,ilt2 
^Qi=K2Hi+Ksi + 


.= E 
^0,-^1,2^^ 


^ /k ^KAg + og H-——, K^Xi + Kti + riij . ( 161 ) 

Note that the s and ng dependence of the right hand 
side comes only from og. Also, after converting the sum 
J2no Sao’ ^be Ri dependence is gone. So, we write 
^Q“=K4*-+Kii,+r. xaonin 2 .iit 2 ^ Observe that 

^Q°’^K 2 A+Ks +r to depend on nine Zx-valued 

parameters, ngp^ 2 , ^ 1 . 2 , si, 2 , ?’i .2 After the reorganization 
we have just done, we lost Si^ 2 ) and we now only have 
six Zk parameters, ni^ 2 ) ^i, 2 Xi ,2 and og. While og is not 
Zx-valued, we can shift og such that 

( 162 ) 


where the second term takes values 0,... K — 1. Then, 

vlao]nin2,tit2 _ \ ' 

rir 2 / , 

Aq,Ni^ 2 ^^ 

X /k ^KAg + [og] H-—, K^Xi + Kti + Ui^ . 

(163) 

Observing that depends on [ag],ni_ 2 Ai ,2 

only through ngp _2 defined in Eq. (163), rewriting Eq. 
(163) in terms of hg 1 2 completes the derivation of Eq. 
(153). 


3. Modular transformations 

We now discuss the modular properties of the partition 
functions. Under the U 2 transformation, the zero-mode 
partition functions are transformed according to 


(172X)"°’’'^"^ = g“ j^no"in2 


(164) 


On the other hand, under the transformation, the 
partition functions are transformed as 


■/ Y\nQnin 2 


{U[X) 


1 


g|(¥("on'i-l-ni"o)x"o"l"2 

,1 

(165) 


where fig = Kng -|- r 5 ^n 2 — 

Observe that, upon the U[ transformation, partition 
functions with new parameters r( and generated. 
Since ri and r 2 are the given quantum numbers from 
the 0 term, the action of modular transformations is not 
closed. 


V. COUPLING TWO BF THEORIES - 
THREE-LOOP BRAIDING STATISTIGS 

In the twist (124), the integers Qi ^2 are fixed and 
treated as a background. I.e., 0 is a non-dynamical field. 
We have seen that the surface partition functions do not 
form a complete basis under modular transformations. 
To circumvent this issue, one may consider to treat Qi 
and <52 as dynamical variables, which may come from 
another copy of the BF theory. In this section, we will 
discuss two copies of the BF surface theories which are 
coupled via cubic terms. 

Let us start from two decoupled copies of the BF sur¬ 
face theories. Let Mg_i _2 and Qg,i ,2 label different twisted 
sectors of the first and second copy, respectively. We con¬ 
sider to twist these quantum numbers by 


Mg ^ Mg -h 

Qo Qo + 


Q X M 
K ’ 
M xQ 


ttg = K X (integer) -|- [ug] 


K 


(166) 
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Here, unlike Eq. (124), both M and Q are dynamical 
variables. 

In the next section, we start by introducing an (3+l)d 
bulk field theory, Eq. (167), or its alternative form (173), 
which realizes precisely the twist (166). We will analyze 
the modular properties of the resulting zero mode parti¬ 
tion functions at the surface. The oscillator part of the 
partition function is simply given by the partition func¬ 
tion of the two decoupled copies of free boson theories. 
By computing the S and T matrices acting on the zero 
mode partitions, we argue that the action (167) realizes 
three loop braiding statistics. While we were finalizing 
the draft, a preprint 54 appeared where the similar bulk 
actions were discussed and conjectured to realize three- 
loop braiding statistics. 


A. The bulk field theory 


a. The cubic theory Let us motivate the twist (166). 
We propose to work with the following bulk action: 


Shulk — 



K 

27r 


^ijb^ A da"^ 


+ A A da^ + A A da^ 

47r2 47^2 


— 5i.jb^ A Jqy — Sija^ A J^p , 


(167) 


where I,J = 1,2 and pi _2 are, as the level K, constant 
parameters of the theory. Similar action has been dis¬ 
cussed in Ref. 55-57. The equations of motion are 




K 


Pi 2 


——db a doj 

27T 4:7T^ 

-^a^Ada^ + p.da^Aa^ = Jl, 

27r2 47r2 

—db^ + A da^ 

2 tt Att^ 

- A da^ -I- -p^da^ A . 

27r2 47r2 


(168) 


As in our previous discussion in the BE theory with 
and without the theta term, let us consider a fixed, static 
quasiparticle and quasivortex configuration and integrate 
the equation of motion over space. By solving the first 
equation of motion as = (27r/K)(d“^ J^^), plugging the 
solution to the second and the third equations of motion. 


and integrating over space, 

+ ^ A Jgp, 

+ ^ A) A + Ap’ (169) 

where note that in the static configurations considered 
here, Jqv is a delta function one form supporting a spatial 
loop, whereas Jqp is a delta function three form support¬ 
ing a spatial point. Correspondingly, d~^Jqv is a delta 
function 0 form supporting a three dimensional manifold. 
The contributions to db^ is coming from quasivortex 
loops, J^{d~^ J^y)AJqj,, are given in terms of their linking 
number. 

Considering now the specific geometry E = x 
with the boundary (surface) dH = T^, we can derive 
the quantization rule of the zero modes of the bound¬ 
ary fields. Using the Gauss law constraint to write the 
boundary conditions in terms of and the bulk equa¬ 
tions of motion translate in to 



With p;^ = p 2 = K, these correspond precisely to the 
twist (166). 

Note that if we naively gauge transform as b^ ^ b^ + 
dC^ and ^ a + d(p\ we find that the theory is not 
gauge invariant. Moreover these gauge transformations 
are not generated by Gauss constraints. We propose the 
following alternative gauge transformations: 

b^ -A b'^ = b^ + d(^ — (a^ A d(p^ + dp'^ A a^) , 

5^ — b'^ = b^ + d(^ — (a^ A dp'^ + dp^ A a^) , 

-A a'^ = -\- dp^. (171) 

Therefore the action with cubic terms in Eq. (167) is 
gauge invariant. On the other hand, as for the coupling 
to the sources, by demanding the gauge invariance, we 
can read off the conserved currents, which are modified 
due to the presence of the cubic terms and the modified 
gauge transformations. 

On an open manifold, the action picks up a gauge 
anomaly on the boundary under these gauge transfer- 






19 


mations 


where 


^bulk\P 1 ^ 

+ ^ [ SijdC^ A a'’ 
JdM 


+ A dip^ A + p2(i(A^ A dtp^ A a^). 

JdM 

(172) 


This anomaly then must be compensated by an appro¬ 
priate boundary field theory. 

b. The alternative quadratic theory Instead of tack¬ 
ling the cubic theory (167) and the corresponding surface 
theory, as in our discussion in the BF theory with theta 
term, we consider an alternative form of the theory. We 
note that the equations of motion (168) can be derived 
from the following alternative action: 


<?' = 
^bulk 


^ J dijb^ A da"^ 


P2 1 
-.a 


27rK 


62 + 

27rK 



A 

A 


5ija'^ A Jqj 
A Jlv 
A 


(173) 


Unlike Sbuik, this theory is quadratic. Integrating over 
and , one obtains the effective action of the currents 


J I?[a^ (174) 

where 

5e// = -^ 

+ Pi J (d-^Jqv) A id~^Jqv) A Jg^ 

+ (^) P2 J{d-^Jqv)/\id-^Jlv)/\Jly (175) 

The first term in the effective action describes, as in the 
ordinary BF theory, the quasparticle-quasivortex braid¬ 
ing statistics while the second and third terms include 
interactions among three quasivortex lines. 

From the coupling to the currents, we read off the Wil¬ 
son loop and Wilson surface operators in the theory: 


■= A a^, (178) 

27rK 

and as before the repeated capital Roman indices are not 
summed over. Note also the triple commutator among 
fg is computed as 


fgB^Jg,B^\jg„B^ 


K3 


I{S#S',S"). (179) 


To make a comparison between the cubic and quadratic 
theories, in the cubic theory, the canonical commutation 
relations differ from the ordinary BF theory, while they 
remain the same in the quadratic theory. In fact, in the 
cubic theory, the commutator among fields generates an¬ 
other field, [6, b] ^ a, schematically. On the other hand, 
the set of Wilson loop and surface operators in the cubic 
theory is conventional (i.e., identical to the ordinary BF 
theory) while it is modified in the quadratic theory as in 
(176). In spite of these differences, the algebra of Wilson 
loop and surface operators of the two theories appear to 
be identical. Therefore, we argue that the two theories 
are equivalent. In the following, we will proceed with the 
quadratic theory. 

c. the quantization rule of the zero modes We now 
derive the compactification condition of the boundary 
fields from Eq.(176). In the presence of a boundary and 
using = df^ and = dp^, the surface operators re¬ 
duce to 


expiTO 



Pj 

27rK 


dp^ A dp^ 


=exp im 



Pi 

27rK 


p^ A dp^ 


(180) 


where the boundary of the world sheet is on the surface. 
We now consider the case where dS is along the x- or y- 
cycles on the surface. Recalling the mode expansion 


= «o 

c;(r) = 




Ih 

+ 


+ 


Pjy 

i?2 

% 


271 Rj 27ri?ii?2 


x5j^2 + ■ 


(181) 


exp 

im / 

, exp in 


J L . 

Js ^ 


A ad 

27rK 


, (176) 


where n and m are integers, L and S are arbitrary closed 
loop and surfaces, respectively, and we introduced the 
notation 1 = 2 and 2 = 1, and the repeated capital Ro¬ 
man indices are not summer over here. These operators 
(or rather their exponents) satisfy 




27ri 

i<r 




2ipj 


d Is#s' + 


J , 2ipj^/j 


K2 


Is#S’ 

(177) 


the zero modes enter into the integral /gg[C^ + 
{pj/27rK)p^ A dp^] through the combinations 

(182) 

We thus conclude 






(183) 


are angular variables, where 

u^° = (1,0,0)^, 

= (Pr/3f/K, 1,0)^, ^ 0,1)^. (184) 
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The rest of the discussion is essentially identical to the 
analysis made in Sec. IV B 1. We recall the commutation 
relations among zero modes 

K>/5o] = (185) 

the following linear combinations of the zero modes are 
integer-valued 


Under the U 2 transformation, the product 
^norain 2 ^rorir 2 jg invariant up to a phase, 

_ ^-^('noni+fori)-^(rin2-r2.nx)(ni-ri)j^nonin2j^rof-if2 

(193)’ 

where we have introduced 


KwI^PI = ml, ml e Z, (186) 

where 

wl = (l,-prl3i/K,-pTPl/K^ , 

«;( = (0,1,0), a;2" = (0,0,1). (187) 

Inverting this relation, the eigenvalues are given by 

KPI = Mo - ^(Q X M), 

K/3i=M2, Kpl=-Mu 

K/3o'= Oo - ^(M X Q), 

K /32 = Q 2 , KPl = -Qi, (188) 

where M^ and are integers. 


B. The surface partition functions 

With the twist (166), the two copies of the surface 
theories are coupled together. The partition functions 
are given by 


ynonin2 yrorir2 
^rir 2 ‘^nin 2 


(189) 


ho = Pq — r 2 ti -b r 2 Si mod K, 
f 0 = lo -b 77.2^1 — n 2 Si mod K. (194) 

For the above equation, if we write down the phase in 
terms of uq and Tq, it will be independent of ti and s^. 
In other words, for two different if they 

have the same ho, ho, n-i, ri, 712 and r 2 , the phases they 
acquire under the U 2 transformation are the same. This 
motivates us to combine these partition functions and 
define, for fixed ho, ho G Zk, 


Y- 

Ar 


E 


j^nonin 2 j^rorir 2 


(195) 


where the sum is taken over all quartets (ti, 2 )Si^ 2 ) giv¬ 
ing rise to given ho, ho. Observe that is labeled 

by two Zk X Zx-valued quantum numbers, and one Zk- 
valued quantum number. On the other hand, 
depends on six Zk x Zx-valued indices. There are K® 
sectors. From Eq. (193), it is straightforward to read off 
the transformation of y under the U 2 transformation: 

_ ^-^(noni+fori)-j^{rin2-r2ni){ni-ri)^nonin2 ( 2 . 96 ) 


where, as before, we decompose the quantum numbers as 

M^ = K7V^-bn^, = 0,1,... ,K - 1, 

Q^=KRf 2 +r^, = 0,1,... ,K - 1, G Z, 

(190) 

and noted, following the discussion in Sec. IV B 2, the 
partition functions depend only on the fractional parts 
of M^/K and Q^/K. Following Sec. IV B 2 further, we 
can write the partition functions as 


As for the U[ transformation, the product 
X^on^^ 2 xfofif 2 transforms under U[ as 


where the phase 9 is given by 


.|2K4 


(197) 


i.T-bieZK 


ynonin2 yrorir‘2 
^rir 2 ^n\n 2 


_ ^ ^ j^nonin2 j^rorir2 

i 1,2 ,-S1,2 ^ 


(191) 


where 


tiq := no + s X n + r X t mod K, 
Tq := ro -b t X r -b 71 X s mod K, 


Tig = Kpq -b (r X 7l) 
hi = Kti -b Til 
h2 = Kt2 + 712 


ho = Krg -b (ti X r) 
hi = Ksi -b ri 
h2 = Ks 2 -b 72 


( 192 ) 


on’on'R'oK ^ ^TTTioTi'i 27r7ii7i() 27rror'i 27rrir() 

K2 K2 K2 K2 

Stt 

= ^{non[ + hgTii -b foY + hgri) 

O-TT 

+ J^[(7 X 7l)(7l'i -r'l) 

+ {r[n 2 -r 2 n[){ni-ri)]. (198) 

To derive this result, we note that 712 and 72 are invariant 
under the U[ transformation. As in our previous discus¬ 
sion on the U 2 transformation, it is crucial to observe that 
the phase is independent of ti, 2 ,si_ 2 - We are 
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(n'l, n; 712 ) 

(ni,ri;r2) 

(7-(,ri;n2) 

(n'i,ni;r2) 

27rnon\ r-| 

2'Kron\ r-| 

ATTnor'-t r-| 

ATTron\ m 






(n'l, n; 712 ) 

(n'l, n; 72) 

(r(,ri;n2) 

(n'l, m; 72 ) 

27rq2n2n'j^ri 

27rqi r2 n'-^ ri 

47rqin2'Tj'ri 

47rq2r2n^ ni 

K2 

K2 

K2 

K2 


TABLE I. The braiding statistical phases (the second line) 
for the braiding processes between loop a and loop b with 
base loop c linking both of them (denoted by (a, b; c) in the 
first line). Here, a, b and c are the quantum numbers for loop 
excitations. 


thus led to consider the partition functions Xf°ri^ 2 ^ 
fined in Eq. (195), which transform, under the U[ trans¬ 
formation, as 


(TJ/ ^\'nonin2 

1 A./rQrir‘2 


|t|2K2 


”o."l,2.ro.ri,2 


Y y '5 


(199) 


Summarizing, the modular S and T matrices are given 
by 


<Sn 


— S , p-^in'o^2-non[+far2-for[) 

X [("i+ri)("2r(+n(r2)-2n2n(ri-2nir2r(] 


Tn, 


= 

2 

X e 


-{noni+fori)+^{rin2-r2ni)(ni-ri) 


( 200 ) 


where = (ho,ni,n 2 ). Observe that a quick way to 
obtain this three loop braiding phase is to replace uq —>■ 
no -b (r X n)/K, ro —> tq -b (n x r)/K, in the S and T 
matrices for the surface of the BE theory (Eq. (113)). 

The first exponential in the S matrix, 
^- 2 Trz{ngn 2 )/K^ |.]^g gj.g^ term in Eq. (198) 

represents the particle-loop braiding phase, which 
exists also in the ordinary BE theory. On the 
other hand, the second exponential in the S ma¬ 
trix, -l/K", and the second 

term in Eq. (198) describes a topological invariant 
which can be considered as the higher dimensional 
generalization of the linking number of closed lines (in 
three dimensions), and is also related with the three-loop 
braiding process.More precisely, from the second 
term in Eq. (198), one can extract three-loop braid¬ 
ing statistical phases. Eor example, the phase factor 
g27ririn2ni/K included in Eq. (198) can be interpreted 
as the three-loop braiding statistical phases associated 
to two loops running in the a;-direction with quantum 
numbers ri and n'l in the presence of a base loop 
running in the y-direction with quantum number n 2 
(Table 1). The three-loop braiding statistics encoded 
in the 5-matrix can be further understood through 
dimensional reduction discussed below. 

As for the T matrix, the first phase factor 
^ 2 m(noni+rori)/K jg pj-Qposed to be the topological spin 
for the composite particle-loop excitations in the BE 


TABLE 11. Same as Table 1, but for generic values of the 
parameter qi ,2 = 0,1, 2,..., K — 1. 


theory. On the other hand, the second phase factor 
g27ri(rin2-i-2ni)(rai-ri)/K^ Considered as the topo¬ 

logical spin for the loop excitations with a base loop 
threading through it. For instance, e 2 ’r*rira 2 ni/K j.gp_ 
resents the topological spin for the loop excitation with 
quantum number (ri, ni) threaded by the loop excitation 
carrying quantum number 712 - 

These results extracted from the boundary S and T 
matrices, (200), are consistent with the previous bulk 
calculations in the literature.In particular, in Ref. 
34 the S and T matrices in the bulk are calculated in 
the basis that is constructed from the so-called minimum 
entropy states (MESs) on the bulk spatial three torus. In 
Ref. 36, the bulk S and T matrices were constructed for 
Zatj X ZjVj X gauge theories. 

Several comments are in order: 

(i) The entropy St computed from these characters 
and the modular S matrix shows, in the limit Ri/Rq —>■ 
00 and R 1 /R 2 —>■ 00, the asymptotic behavior St = 
—21nK -b • • •, where • • • is the term proportional to the 
area of the surface. I.e., the constant piece in the (entan¬ 
glement) entropy is the same as the two decoupled copies 
of the BE theories. 

(ii) For (3 + l)d topological phases (gauge theories) 
with Zk X Zk gauge symmetry, we expect there are 
(at least) different topological phases that are differ¬ 
entiated by their three-loop braiding statistical phases. 
This is expected from the group cohomology classification 
(construction) of SPT phases; from H^[Zk x Zk,U{1)] = 
Zk X Zk, we expect there are at least different SPT 
phases in (3-|-l)d protected by unitary on-site symmetry 
G = Zk X Zk. Once the global symmetry is gauged, 
these different SPT phases give rise to different topo¬ 
logically ordered phases which are differentiated by the 
three-loop braiding phases.®®’^® The model we studied in 
this section, the two copies of coupled BE surface theo¬ 
ries, corresponds to the surface theory of one of the 
topological phases. The surface theories of all the other 
topological phases can be obtained by tuning the coeffi¬ 
cient in the coupling terms. In our model, the coefficient 
Pi and p 2 in front of the cubic terms are chosen to be 
K. In general, they can take value qiK and q 2 K with 
91.2 = 0,1,..., K — 1®° which will lead to different 
topological phases with different S and T matrices. The 
three-loop braiding phases will be slightly modified and 
are shown in Table II, which are consistent with Ref. 35. 

Observe also that for G = Zk, H'^[Zk, b^(l)] = 0, i.e., 
there is no non-trivial SPT phase protected by G = Zk 
symmetry. Hence, there is essentially only one topologi¬ 
cally ordered phase with Zk gauge group, whose surface 
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is described by the one-component surface theory studied 
in Sec. III. On the other hand, the two-component surface 
theory studied in this section allows richer possibilities. 

(iii) An insight on the three-loop braiding statistics 
phase can be obtained from dimensional reduction. For 
the trivial two-component BF theory, there is only a non¬ 
trivial particle and loop braiding phase described in Eq. 
(113). This model, after dimensional reduction, reduces 
to the D{Zk X Zk) quantum double model with the K- 
matrix given by Kcrj, 0 Kctj, . 

For the topological phase with non-trivial three-loop 
braiding statistics phase, the dimensional reduction is 
more interesting. Here, we consider the simplest non¬ 
trivial example with K = 2. We perform dimensional re¬ 
duction on U[ and T defined in Eq. (199) and Eq. (200). 
When we do so, we need to fix the quantum numbers 
n 2 and r 2 - For example, if we take n 2 = 0 and r 2 = 0, 
i.e., there is no third loop connecting the first and second 
loops, the S and T matrices after dimensional reduction 
are the same as those for the two copies of the toric code 
model. 

On the other hand, if we take n 2 = 1 and r 2 = 1, the 
dimensional reduction results in the S and T matrices 
given by 


^rii ,ri 

q-n'-y- 

/ rii ,ri 


— ^„Triinon[+ngni+ror[+r'gri)+Tri(ri-ni)(n[-r[) 

~ 4 

= ( 201 ) 


This indicates that the (2+l)d topological order de¬ 
scribed by the K-matrix 


/ 2 2 -2 0 \ 
2 0 0 0 
-2 0 2 2 
V 0 0 2 0/ 


( 202 ) 


By an 5'L(2,Z) similarity transformation, this K-matrix 
is equivalent to K = 2az © 2crz, which represents two 
copies of the double semion model. 

Similarly, if we choose (n 2 ,r 2 ) = (1,0) and (0,1), the 
corresponding (2 + l)d topological order is described by 
the K-matrix 


and 


/ 0 2 -1 0 \ 
2 0 0 0 
-10 2 2 
V 0 0 2 0/ 


/ 2 2 -1 0 \ 
2 0 0 0 
-10 0 2 
V 0 0 2 0/ 


(203) 


(204) 


respectively. To summarize, after dimensional reduction, 
the original (3 + l)d topological order with non-trivial 


three-loop braiding statistics “splits” into four differ¬ 
ent (2 0 l)d topological order, which are controlled by 
the quantum numbers n 2 and r 2 . This result seems to 
be related with the group cohomology classification of 
symmetry-protected topological (SPT) phases in (2 + l)d 
with G = Z 2 X Z 2 symmetry, i.e., H^[Z 2 x Z2,17(l)] = 

Z 2 X Z 2 . 


VI. DISCUSSION 


Let us summarize our main results. 

- In the (30l)d BF theory, we have demonstrated, 
through explicit calculations in the boundary field theo¬ 
ries and by comparisons with known bulk results, there 
is a bulk-boundary correspondence in (30l)d topological 
phases. In particular the modular S and T matrices are 
calculated from the gapless boundary field theory and 
shown to match with the bulk results. 

- The surface theory of the (30l)d BF theory with 
the theta term is introduced and solved. The action of 
the modular S and T transformations on the partition 
functions is calculated. It is shown that the partition 
functions do not form the complete basis under the mod¬ 
ular S and T transformations. 

- Finally, we propose a (3+l)d bulk field theory 
with cubic coupling that may realize three-loop braid¬ 
ing statistics. We discuss the twist that the cubic term 
of the field theory adds to the zero modes. By consider¬ 
ing the alternative form of the bulk and boundary field 
theories, in which the quantization rule of the zero modes 
is twisted, we computed the surface partition functions, 
and the S and T matrices are constructed. 

These results extend the well-established bulk¬ 
boundary correspondence in (2+l)d topological phases 
and their (l+l)d edge theories. Our approach from the 
surface field theories provide an alternative point of view 
to (3+l)d topological phases, and to recently discussed, 
novel braiding properties, such as three-loop braiding 
statistics. 

There are, however, still some aspects in the (20l)d- 
(l0l)d correspondence, which we do not know if have an 
analogue in the (3+l)d-(2+l)d correspondence. For ex¬ 
ample, in the case of the bulk-boundary correspondence 
connecting (2+l)d topological phases and (l+l)d edge 
theories, that the edge theories are invariant under an 
infinite-dimensional algebra seems to play a significant 
role: the Virasoro algebra or an extended chiral algebra 
of (l+l)d CFTs faithfully mirrors bulk topological prop¬ 
erties of (20l)d bulk phases. On the practical side, that 
edge theories enjoy an infinite-dimensional symmetry al¬ 
gebra provides many non-trivial solvable examples. For 
our example of (2+l)d surface theories of (3+l)d topo¬ 
logical phases, on the other hand, we did not make use 
of such infinite-dimensional symmetry. In fact, the sur¬ 
face theories studied in this paper are not conformal field 
theories. For example, the two-point correlation function 
of the boson field r)(j){t', r')) in the free boson the- 
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ory in (2+l)d decays algebraically. This should be con¬ 
trasted with the logarithmic decay of the corresponding 
correlator in the (l-l-l)d compactihed boson theory. As a 
consequence, the correlation functions of the bosonic ex¬ 
ponents exp[im(j){t, r)] (m € Z) do not decay algebraically 
in the (3-|-l)d free boson theory. Whether or not there 
exists a unified held theory framework in (2-|-l)d held 
theories that strongly resonates with topological proper¬ 
ties of (3-|-l)d bulk topological phases requires further 
investigations. 
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Appendix A: The surface theory of the BE theory 
with theta term 


In this appendix, we go through canonical analysis of 
the surface theory of the BF theory with theta term. It 
is described by the Lagrangian density 


C = 


2 ^ Ci j di Q + ^^2 ^ij di&dj if 


dt<f 


I 

2Ai 


(etjd.Cjf - (AI) 

Za2 


In this theory, physical observables are bosonic exponents 
exp[im(p{t,r)], and Wilson loops expim dx^Q{t,r), 
where m € Z, and C is a closed loop. The boson Held (p 
and gauge fields Q are compactified accordingly. 

The canonical commutators are 


[cijdiCjit, r), eimdiCmit, r')] = - r'), 

O'jrj 

[ip{t, r')] = (A2) 


The mode expansions consistent with the equations of 


motion are given by (only the oscillator parts are shown) 


p(t, r) = 


1 


1 


,2k-i 




V2K2AiRii?2 yA(^ 

jt(k)e*A+(k)* + a(_k)e*^-«t 

A) ^ 1 


E 




87r2i?ii?2 VA(k) 

tA+(k)a^(k)e*^+«* +zA_(k)a(-k)e*^-(^)‘' . 

(A3) 


where 


S(k) = 


K2AiA 2^ 
A(k)2 = S(k)2 -k A(k)2 

I / 


K2Ai 


A2 


K2Ai 

A±(k) = -S(k)+A(k) 
where 


) hkj 


K2Ar 


jjij _ I —Q1Q2/R1R2 

\—Q1Q2/A1R2 QilRi , 


(A4) 

(A5) 


We now attempt to derive the zero mode quantization 
rule solely from the surface theory. The canonical com- 
muators among the zero modes are 


-1 




(A6) 


where 



( K 

0 

0 ^ 


( k 

0 

0 \ 

M = 

-VQ 2 

0 

TK 

, M-i = 

pQi 

K2 

0 



^+p(5i 

-K 

0 J 


PQ2 
\ K2 

1 

K 

0 J 


(A7) 

The commutators take a canonical form in the rotated 
basis 


i3x = l3„M^x ^ [af^,i3x\ = iS^x- (AS) 

Thus, from the compactification condition on a^, the 
eigenvalues of j3 are integers: 

= wa e Z. (A9) 

This means the eigenvalues of /3 are 

Renaming integers by mo —>■ no, mi —> —n2, m2 —>■ n-i, 

(/3o,/3i,/32) = (^-j^(Qxn),^,^). (All) 

This is consistent with the quantization rules derived 
from the equations of motion in the bulk field theory. 
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Appendix B: 5-function forms 

In this appendix, we summarize the properties of 6 - 
function forms. For an n-dimensional submanifold Af of 
a U-dimensional maniofold A4, we define a {D — n)-form 
So-niAf) by 

f An = f So-niAf) A An, ^An, (Bl) 

JM 

where is an arbitrary n-form on AI. If we flip the 
orientation of N, 

So-ni—Af) = —5D-n{Af)- (B2) 

More generally, for oriented submanifolds Ni, 

=Y,c^5{M,) (B3) 

i i 

where is a coefficient. 

The exterior derivative acts on the delta function form 
as 

5D-n+i{dAf) = {-l)^-^+^d5D-u{Af). (B4) 


Let Ml and be a submanifold of AA with dimensions 
TT-i and 712, respectively. Define d as 

d = ni + n2 — D. (B5) 

When d > 0, Afi and can have a d-dimensional inter¬ 
section within AA. By properly defining an orientation, 
we define the intersection of Afi and A/2, I = Afi#Af 2 - 
The orientation of I is defined to be consistent with 

do-dil) = <5D-ni(A/l) A l5z5-ni(A/2). (B6) 

If the two submanifolds have complementary dimensions, 
n-i + n 2 = D, they intersect at points. Then, the inter¬ 
section number 

I{Afi,Af2) = J SniAfi) A fe-„(Af2) G z (B7) 

counts the number of intersection points. 

The linking number of two submanifolds, C and Af, is 
defined when dim£ -|- dimAf -I- 1 = D. By considering 
an auxiliary manifold satisfying dAfi = C, the linking 
number is given by 

Lk{C,N) = I[Ni,N). (B8) 
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